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ON FLOW OF ELECTRICALLY CONDUCTING FLUIDS OYER 
A FLAT PLATE IN THE PRESENCE OF A 
TRANSVERSE MAGNETIC FIELD 1 

By Vehnon J. Rossow 


SUMMARY 

The use of a magnetic field to control the motion of electrically 
conducting fluids is studied. The incompressible boundary- 
layer solutions are found for flow over a flat plate when the 
magnetic field is fixed relative to the plate or to the fluid. The 
equations are integrated numerically for the effect of the trans- 
verse magnetic field on the velocity and temperature profiles, 
and hence, the skin friction and rate of heat transfer. 

It is concluded that the skin friction and the heat-transfer 
rate are reduced when the transverse magnetic field is fixed 
relative to the plate and increased when fixed relative to the fluid. 
The total drag is increased in all of the cases studied. 

INTRODUCTION 

It has been said that a fluid at a very high temperature is 
like a universal solvent which cannot be contained. A 
possible method of containing this fluid is suggested when 
it is noted that at such a high temperature it would surely 
contain ions and quite probably also free electrons. The 
fluid would then be an electrical conductor. The invisible 
hand of electrical and magnetic fields can then be used to 
induce forces on the fluid such that it is prevented from 
coming in direct contact with a wall which it would dissolve. 

A somewhat similar technique has been in use for some 
time in the metal purification industry. It employs a high- 
frequency magnetic field which causes eddy currents in a 
lump of molten metal which in turn react with the imposed 
magnetic field. The metal is thereby suspended in space 
if the imposed magnetic field is made strong enough. 

Another example is the so-called “perhapsatron” described 
briefly in reference 1. A gas in a doughnut shaped container 
is heated to a high temperature by an electrical current dis- 
cliarge. Through the use of the interaction of the resulting 
ion current and a magnetic field, the hot gas is prevented from 
coming in contact with the surface of the vessel. An appli- 
cation of similar principles was used in thermonuclear fusion 
experiments in the Soviet Union. The techniques and 
results described very briefly in reference 2 indicate that it 
was possible to keep the hot fluid from the walls and to 
concentrate the hot gases quickly so as to generate a focusing 
shock wave. 


These examples indicate the possibilities which may be 
realized through the use of a magnetic-field force on a flowing 
conducting fluid. A number of situations exist in aero- 
dynamics wherein a magnetic or electrical field might be 
used to alleviate high convective heat-transfer rates to a wall. 
Such problems arise in the flow of air in the boundary layer 
and in the vicinity of stagnation regions of an aircraft moving 
at very high speeds. If the velocity is high enough the air is 
ionized and, therefore, electrically conducting. See, for 
example, references 3 and 4. Other examples are those 
associated with the flow of the combustion products (which 
are hot and generally electrically conducting, refs. 5 through 
7) in the propulsion units of aircraft. 

The extent to which a given heating problem can be 
alleviated is as yet difficult to determine. It would be well 
if a few theoretical results were available to evaluate the 
effect of a magnetic field on the drag and heat-transfer rate. 
An attempt to include all the aerodynamic features of the 
flow of air over or inside an aircraft plus the magneto- 
hydrodynamic (or more correctly, magneto-aerodynamic) 
effects would render the problem so complicated as to make 
its solution difficult if not impossible. It is felt that a few 
simple basic theoretical solutions would point out the 
advantages or disadvantages and yield a rough estimate of 
the various quantities entering the problem. It is the 
purpose of this report to present the results of several such 
solutions which, it is hoped, will extend the information 
already available so that such estimates can be made. 
Several applications of the results are indicated. 

PRINCIPAL SYMBOLS 




b,B 

magnetic induction, lines/sq in. 

Cf 

local skin-friction coefficient 

c v 

heat capacity of air at constant pressure, 
Btu/slug °F 

E 

electric field intensity, volts/in. 

E 

/ U^ 

total energy, C p T-\—^ 

erf 

error function 

erfc 

complementary error function, 1-erf 

fo 

dimensionless Blasius boundary -layer 

stream function 


i Supersedes NAOA TN 3971 by Vernon J. Rossow, 1957. 
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dimensionless stream function, equation 
(26) 

force, lb 

local heat-transfer coefficient, Btu/sec ft 2 
°F 

£ 

magnetic intensity, —> (ampere) (turns) /in. 
electric current, amp 

electric current density in the fluid, amp/sq 
in. 

electric current density associated with 
generation of magnetic field, amp/sq in. 
thermal conductivity, Btu/sec ft 2 °F/ft 
characteristic length 

magnetic parameter, per in. 

<rB* 

1 per sec 

P 

magnetic Hartmann number, -^/ ^ 

pressure, Ib/sq in. 

RVZ 

P U 

Reynolds number 
Laplace transform variable 
time, sec 
temperature, °R 

velocity components in x,y,z directions 
Laplace transform of velocity u (see eq. 
(14)) 

velocity, ft/sec 

coordinate axes, x alined with free stream 

excess charge density, coulomb/cu in. 

magnetic viscosity, ^ _ 

coefficient of viscosity, slugs/ft sec 
magnetic permeability, (volt) (second) / 
(ampere) (inch) 

kinematic viscosity, ^ ft 2 /sec 


density of fluid, slugs/cu ft 
conductivity, mhos/in. 
shear stress, lb/sq in. 
stream function 

SUBSCRIPTS 

basic quantity 
free stream 
magnetic induction 
electron 

at plate, j?=y= 0 

SUPERSCRIPTS 


vector quantity 
( )"' derivatives 


HISTORICAL REVIEW OF MAGNETO-HYDRODYNAMICS 

The earliest known published works treating a problem in 
the flow of an electrically conducting fluid through a mag- 
netic field axe those of Hartmann and of Hartmann and 
Lazarus in references 8 and 9. Since that time a number of 
theoretical and experimental studies have been carried 
out. Some of these pertain to the flow of a conducting fluid 
but a larger number relate to the dynamics of ionized clouds 
or stars. No attempt will be made here to include a com- 
plete discussion of the latter group. However, much can 
he learned by studying the papers on the flow of conducting 
fluids. These papers will be divided into several groups 
according to the type of problem treated. 

1. Flow in a channel (refs. 8 through 19). 

2. Flow about bodies of various shapes (refs. 20 through 

22 ). 

3. Related papers on astronomy (refs. 23 through 26). 

The various papers ■will be reviewed briefly on a group basis 
and reference tb particular papers given only in special 
instances. 

FLOW IN CHANNEL 

Following his invention of the electromagnetic pump 2 in 
1918, Hartmann, reference 8, developed equations which 
describe the steady-state flow inside a channel in the presence 
of a magnetic field or a magnetic and electrical field. Inter- 
est was stimulated in problems of this type when it was found 
in reference 9 that a turbulent stream could be stabilized to 
the extent that it could be forced to return to a laminar flow. 
Such a sequence of events is opposite to the usual form of 
transition when a magnetic field is not present. Further 
investigations by other authors, references, 10 through 19, 
substantiated these results and led to the introduction of 

terms such as “magnetic viscosity,” A=— 1 ■ _: > “magnetic 

Hartmann number,” M=Bl x - and 

V a pu 

A rigidity of the fluid is brought about by a coplanar 
magnetic field so that stabilization is achieved because the 
conducting fluid experiences a force which resists the motion 
across the magnetic lines of force (ref. 13). For oxample j 
if the main fluid flow direction is parallel to the magnetic 
lines of force, no interaction takes place unless somo devia 
tion or instability arises. In this event, a restraining force 
will he induced which is proportional and opposite to the 
velocity component perpendicular - to the lines of force, 
tending to damp out the initial instability. 

When the magnetic field is perpendicular to the flow 
direction, a change is brought about in the velocity distribu- 
tion. This may make the fluid flow more stable or unstable 
to transition to turbulence. In the case of flow of an elec- 
trically conducting fluid through a transverse magnetic field 
in a two-dimensional channel (refs. 9 and 14), the change in 
the velocity profile had a very favorable effect on transition. 
Such is not always expected to be the case in other flow 
problems. 

3 Electromagnetic pumps are used in atomic energy plants to transfer radioactive 
solutions from one place to another (e. g., see p. 32 of Atomic Power — a ScicnttGo American 
boot, Simon and Schuster Inc.. New Yort, 1965). 
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These two techniques which increase the stability of fluid 
flowing in a two-dimensional channel were studied theoreti- 
cally in references 13 and 14. The transition boundaries 
are given in graphical form for a range of the magnetic 
parameter. It was observed that a considerable increase in 
laminar run could be achieved at the higher values of the 
magnetic parameter. 

FLOW ABOUT BODIES OF VARIOUS SHAPES 

The papers describing the external flow of a conducting 
fluid about bodies are small in number. The three references 
listed, 20 through 22, describe the flow of a conducting fluid 
around a circular cylinder and a sphere in the presence of a 
magnetic field; the third paper treats the stability of the 
flow between rotating cylinders. Streamlines are shown in 
the first two cases and the stability boundary for the latter 
case. Problems in this category are naturally more difficult 
than the unidirectional channel flow problems discussed in 
the previous section. 


where the term BE results from the electrostatic force on the 
excess charges due to the presence of an imposed electric 

field E. The second term describes the force on the fluid 
due to the interaction of the electric current, j, in the fluid 

■4 

and the magnetic induction B. The differential equation 
can be written in vector notation as 


DU BE 
Dt p 


- jX3 + - gradp=yv s f7 
p P 


( 1 ) 


where p is the fluid density, d is the excess charge density, 
and JET is the magnetic permeability. The substantial deriva- 
tive D/Dt is associated with the use of spatial coordinates. 
A simplified version of Maxwell’s equations appropriate for 
the present problem is 

div 27=0; div J— 0; ]iH=B (2a) 


BELATED PAPERS ON ASTRONOMY 

t 

A need for understanding the motion of electrically con- 
ducting fluids (ionized clouds) has been felt by astronomers 
for some time. The dimensions and quantities considered 
are generally not reproducible in the laboratory or in the 
earth’s atmosphere; however, some of the equations and 
results may be applied to aerodynamic problems. 

In reference 23, Batchelor investigated the possibility of 
spontaneous magnetic fields arising in a conducting medium 
as a result of turbulence. He found that unless 4T£<n»>l 
such fields cannot occur. In all foreseeable problems in 
aerodynamics iirp<rv is many orders of magnitude less than 1, 
ruling out spontaneous magnetic fields in a fluid due to 
turbulence. 

The manner and velocity at which pressure waves are 
propagated in an ionized fluid and other cosmic problems 
are treated in references 24 through 26. The results of 
these studies on wave propagation may be useful in studying 
interference problems on high-velocity aircraft. 

MAGNETO-AERODYNAMIC EQUATIONS 

The equations which describe the flow of an electrically 
conducting fluid in the presence of magnetic and electrical 
fields will now be discussed. The equations are simplified 
to the extent that the flow is assumed to be incompressible 
and to have constant properties (C p , p, k) throughout the 
flow field. In all cases except one, the electrical conductivity 
of the fluid is assumed to be constant. 

MOMENTUM-TRANSPORT EQUATIONS 

The momentum-transport equations for the flow of a 
viscous incompressible fluid consist of a combination of the 
force terms arising from the excess charge density 8, and in- 
duced magnetic effect due to the motion of the conducting 
fluid through magnetic lines of force (see, e. g., refs. 23 and 
27) and the usual Navier-Stokes relations (see, e. g., ref. 28). 
In equation form, these force terms are, 

F=dE+]xB 


curl D—i-irJ (2b) 

A reduced form of the generalized Ohm’s law is, from 
reference 29, 

j—a(E+pUXD) (2c) 

The more complete form given as equation 2-22 in reference 
29 includes additional terms accounting for gravity and pres- 
sure gradient effects. Both of these factors will be assumed 
negligible in this analysis. It will also be assumed that the 
conductivity, <r, does not have directional properties. 

The continuity equation is 

div 17=0 (3) 

Equation (2b) expresses the relation between the magnetic 

^ — > 

intensity H and the total electric current density J. The 

distinction between the current densities j and J lies in their 
location and action. Consider first the situation when the 

fluid is stationary and the electric field intensity E is zero. 
The current density in the fluid is then zero. The magnetic 

~4 ^ -) -4 

field strength H 0 or magnetic induction B 0 =pH 0 is a result 
of an electric current J 0 outside of the fluid. The current 

—4 

density J 0 could be thought of as an electric current in a 
coil or electromagnet outside the flow region producing a 
magnetic field of a given number of lines of force per unit 

■4 *4 

area. The relation between J 0 and H 0 is given by equation 
(2b). The basic magnetic field strength will be designated 

as H 0 . 

When a conducting fluid moves through the magnetic 

—4 

lines of force, B 0 , the positive and negative charges are each 
accelerated in such a way that their average motion gives 

- 4 — 4-4 — 4 — 4—4 

rise to an electric current j=crUXB, where B=B 0 +b. The 

—4 

quantity b is the magnetic induction resulting from the elec- 

~4 —4 

trie current j in the fluid. In this analysis b will be con- 
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sidered as a perturbation on the basic field strength and 

• • • — 4 

negligible in comparison with B„; that is, from equation (2b), 


or 


curl (-B 0 +6) =4irp(J 0 -\-j) —4vpuJ 

■4 —4 

curl B 0 =4.7TtiJ 0 


curl b=4irfij 

—4 —4 

where B^b. 

—4 

If an electric potential E were to be imposed across the 

— 4 —4 

flow field, the current density would be changed by j—aE. 
However, the boundary conditions will be taken so that the 

• ■ — 4 —4 -4 

impressed potential, E, is zero. Therefore, j=<rUXB 0 ; that 
is, the electric current is proportional to the voltage generated 
by the relative motion of the fluid and magnetic field. 

The fluid is assumed to be ionized and thereby an electrical 
conductor. However, within any small but finite volume the 
number of particles with positive and negative charges are 
nearly equal. The total excess charge, 6, in this small but 
finite volume will be taken as equal to zero. It has been 
shown in reference 17 that a very small increase in viscosity 
is realized due to the increased ordered-migration velocities 
of positive and negative ions in the presence of the electric 
field. This effect will be neglected and the electrostatic force 

—4 

term, OE will be assumed to be identically zero. 

— 4—4 

The induced electromagnetic force term, jXB 0 will now 
be related to the local velocity vector. The procedure used 
in references 10, 19, and 22 which results in the magnetic 

— rra 

parameter (ratio of electromagnetic and inertia forces) 

and the magnetic viscosity A— ^ 1 _ will not be used. In- 
stead a procedure will be followed which will lead to the 

<?B 2 l 

magnetic parameter, Q==-^j, used by references 7, 9, 11, 12, 
and 13. 

— 4—4 

A relation between jXB and the velocity is obtained from 
equation (2). It will be assumed that: 

(1) Magnetic field lines are perpendicular to free-stream 
velocity. 

(2) Permeability JL is constant throughout the fluid. 

(3) The assumptions already made are that the excess 

—4 

charge density 9 and imposed electric field intensity E are 
assumed to be zero. 


From vector relations, 


(uxB^ XB 0 =- p \u(b 0 -b'^-B 0 (So-tf)] 

— 4 

However, it has been assumed that B 0 is perpendicular to 

—4 

the free-stream velocity U a (assumption (1)). In the prob- 
lems treated, the local velocity direction differs a negligible 

—4 —4 

amount from the free-stream direction so that B o -U=0 and 
the linearized magneto-hydrodynamic force term becomes 


i 




( 6 ) 


The parameter 3 m is related to the viscous Reynolds num- 
ber and the Hartmann number of references 8 through 12 by 

, <xB 0 H oB 0 H 2 1 M 2 
mL ~ pU„- M pUj~Re 


Equation (1) may then be written as 


DU , crB 2 ft . 1 , ' ft 

p - U-\— p grad (7) 

In the development of equation (7), it was assumed that 
the velocity of the magnetic field is zero. Since only the 
inductive force is being considered, the second term becomes 

— — ( U—U B ) for a magnetic field in uniform translation, 

—4 ^ 

where XJ B is the velocity of the magnetic field. 

When electrical currents flow in a plasma it is well to 
define the conditions at the solid boundaries of the flow 
field. The problems treated are assumed to be two-dimen- 
sional with the electric current flowing from and to infinity 
parallel to a flat plate and perpendicular to the stream 
direction. The configuration may also be thought to con- 
sist of flow in a wide channel with the two side walls as 
conductors connected by a conductor of negligible resistance 
located outside of the influence of the flow field. The flow 
in the middle of such a channel should approximate the 
two-dimensional type of flow being analyzed. Such phe- 
nomena as charge accumulation and boundary effects at the 
walls will be assumed negligible. 


(4) Induced magnetic field b is negligible, linearizing the 
induced magnetic force term. 

The first assumption is not highly restrictive since the per- 
pendicular component of the magnetic field is the only part 
contributing a force which will change the laminar velocity 
profile. The other assumptions are reasonable physical 
approximations. 

If these assumptions are used together with equation (2c), 
the induced magnetic force term becomes 

-- {uxB^XBo (5) 


THERMAL ENERGY TRANSPORT EQUATION 

The differential equation describing the relationship be- 
tween the convection and conduction of thermal energy 
and the work done on an electrically conducting fluid in the 
presence of a magnetic field is found by considering the 


* From electric motor equations, it can be shown that 


aT > «- force 
• Z*x velocity 

and therefore 

W \ / -rr is . / magnetic pressure \ /magnetic Reynolds 
pUco \pUa>*/ K * w \dynamlo pressure/ \ number 
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energy entering and leaving an elementary cubic box. The 
result is given by 


(l+D’+d+SKS+S)'} ® 


The derivation of equation (8), with the exception of the 
second term on the right of the equal sign, can be found, 
for example, in reference 28. The term on the left of the 
equal sign takes account of the enthalpy transported by 
the motion of the fluid or by convection. The first term on 
the right represents heat conducted from one element of 
fluid to the next by molecular motion. The expression, 
j^/c, is the heat added by the electrical current produced 
by the motion of the fluid through magnetic lines of force 
or by an imposed electrical field (see eq. (2c)). The re- 
maining terms arise when work is done on the fluid by 
either pressure or shear forces. The symbol ZtyZ)£ is the 
substantial derivative of the particular quantity and V 2 
denotes the Laplacian. 

If the general expression for the electric current density, 
as given by equation (2c), is inserted into the expression 
for the heat added by electrical means, the result is 

yp / — > ~ > -> 

E-\-TJ XB 

It has been assumed that the electric field E is zero and that 
the magnetic field lines are perpendicular to the free-stream 
direction. To a good approximation in the problems to be 
considered, equation (9) can then be written as 



£=<r!7 s 5 0 2 


( 10 ) 


where IP must be written as ( U — U B Y when relative motion 
exists between the magnetic field and the plate. Equation 
(8) is then 


^-m-^+«cw+§f+,{2[(g)’+(|)'+(g)’]+ 

(l+SHI+SKS+S)! <"> 


viscosity, and has the same electrical conductivity through- 
out. A third case, in which the magnetic field is fixed on 
the plate and a compensating pressure gradient exists in the 
fluid, will be shown to be equivalent to the case wherein the- 
magnetic field is fixed relative to the fluid. The problem 
wherein the fluid conductivity varies will not be considered. 
The fluid and plate will be assumed to be initially at rest. 
At time t= 0, the plate will move impulsively with a given 
velocity. 

MAGNETIC FIELD FIXED RELATIVE TO THE PLATE 

At time £<0, the fluid, plate, and magnetic field are as- 
sumed to be everywhere stationary. At time £=0 and for 
all later times the plate and magnetic field are moving at 
velocity u=u a (fig. 1). The problem is to find the velocit 3 r - 
time history of the fluid. 



Figure 1. — Flat plate moving impulsively through a fluid in the 
presence of a transverse magnetic field. 


The velocity in a given y=constant plane does not change 
with x. It can then be reasoned that the vertical, v, and 
transverse, w, velocities are zero or negligible together with 
the pressure gradients in all directions. Equation (7) then 
reduces to 


bu . uB 0 i 


•u=v 


bht 


( 12 ) 


IMPULSIVE MOTION OF A FLAT PLATE 

The impulsive motion of an infinite flat plate 4 in a viscous 
fluid serves as a model for the boundary layer on a semi- 
infinite flat plate. The advantages of working the problem 
of the impulsive motion of a flat plate are that it is simple 
enough to yield a result in closed form and suggests the 
proper choice of parameters to be used on the more compli- 
cated problem of a semi-infinite flat plate. 

The velocity profiles will be found for two cases. In the 
first it will be assumed that the magnetic lines of force 
are fixed relative to the plate, and the second relative to the 
fluid. In both cases the fluid is of uniform density and 


Since the magnetic field is moving and the fluid is initially 
at rest the relative motion must be accounted for. Equation 
(12) is then 


bu . crB 0 i 
bt 1 p 


(u—uj=v 


bht 

&2/ 2 


(13) 


The boundary conditions are 

u=u a , y= 0, f>0 
u=0, y~^> 0, £=0 

The Laplace transform of the velocity u is defined as 



i This problem Is often referred to as the “Rayleigh problem” when the Hold Is Incom- 
pressible and nonconducting. 


e~“udt 


(14) 
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Applying the Laplace transformation (see, e. g., refs. 
30 and 31) to the first term in equation (13) gives 

5? *=+»»■■' I >1" 

= 8U 


If the other terms of equation (13) are treated similarly, 
the transformed equation is 


or 


- (Fu . mni„ 

su=p -j—H • 

ay 2 s 

<Fu -/s+m A 
dy>~ u { v ) 


m x u 


niiu^ 

V8 


(15) 


The solution to equation (15) is the sum of the solution of the 
homogeneous equation plus the particular solution, that is, 


-jn^_ 

s(8+m0 


mjii „ 
«(«+?«!) 



The constant Ci is chosen equal to zero to fit the boundary 
condition that u must be finite at y= <» . The integration 
constant Ci(s) is found from the boundary condition at 
y= 0; that is, on the plate 


The Laplace transform of is ujs. Therefore 


2.4 


1.6 

y 

.8 


0 



(a) Magnetic field fixed relative (b) Magnetio field fixed relative 
to plate. to fluid. 

Figure 2. — Boundary-layer velocity profiles on a flat plate set into 
motion impulsively. 

The skin-friction coefficient is expressed as 


p(&tt J %)Uo 

puJ/2 


so that, from equation (17), 



(18) 


(19a) 


Values of c f are shown in graphical form in figure 3. When 
■rrii is zero, equation (19a) reduces to the skin friction for the 
nonelectromagnetic Rayleigh problem. 


and 



u— 


mju , „ 
^(s-I-TOl) 


-Jba-e 

8+m x 



(16) 


Equation (16) can he inverted with the aid of reference 30. 
Combining several terms gives 


u=u a (l -e-»*erf^=) 


(17) 


The symbol erj denotes the error function of the argument. 
The error function is discussed in reference 31 and tabulated 
extensively in reference 32. Typical velocity profiles are 
shown in figure 2 (a). If m x is set equal to zero, equation (17) 
reduces to the result for the Rayleigh problem; that is, 



u=u„(l-erf^j=) 

The velocity is a function of y/2-fvi only and therefore a 
single profile suffices. However, when a magnetic field is 
acting on the fluid, the velocity profiles are not similar 
because they change with time according to e~ mJ . The fluid 
at infini ty is accelerated by the magnetic field so that the 
entire mass of fluid is accelerated by the magnetic field. It 
is, in fact, accelerated more rapidly than when the'only force 
is the viscous action between layers. 


Figure 3. — Skin friction parameter as a funotion of a magnotlo 

parameter. 


The change in the velocity profile is a result of the force 
exerted by the magnetic field on the fluid. The reaction on 
the unit generating the magnetic field is expressed as 




.F/unit area= aB^iu^—ii) dy=<rB 0 2 e~ mU erf 


1 e 


V 


2 


dy 


(19b) 


The force per unit area goes to infinity as the upper limit of 
integration goes to infinity because the magnetic field extends 
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undiminished an infinite distance from the plate. Thereby, 
an infinite amount of fluid is accelerated, resulting in an 
infinite force. 


MAGNETIC FIELD FIXED RELATIVE TO THE FLUID 

At all times less than zero the fluid, magnetic field, and 
plate are assumed to be at rest. At time f=0 the plate 
begins moving with velocity u=u a but the magnetic field 
remains at rest. The differential equation is the same as 
equation (12) because there is no relative motion between 
the fluid at y — <*> and the magnetic field. 

int . ZPu /nn , 

_ +mi1i== ,_ (20) 

The boundary conditions are 

u=u a aty—0t>0 
u=0 at 0 f=0 

Applying the Laplace transformation to equation (20), 
as was done in the previous section, yields the transformed 
equation 

The solution to this ordinary differential equation is 


u—Ci(s)e 



+ry[^ 


where the constant Oj is set equal to zero because of the 
requirement of a finite velocity at ?/=<». The integration 
constant, Ci(s), is found from the boundary condition at 
y= 0, f>0, as 



Equation (20) can be written as 

_ u - V»+">i -7= 

U=Z 1~L — \ — —e 

(8+mi)— m x 


This equation can be inverted with the aid of reference 30 as 



( 21 ) 


The symbol erfe denotes the complementary error function 
defined as l— erf (see, e. g., p. 370 of ref. 31). The velocity 
is once again dependent on more than one parameter so that 
a single similar profile cannot he drawn. Several profiles 
are shown in figure 2 (b). The fluid at y= is not disturbed 
in this case because the velocity across the magnetic lines 
is zero. In the vicinity of the plate the induced magnetic 
force counteracts the acceleration force of viscosity resulting 
in an increased rate of shear at the walls as expressed by the 
skin-friction coefficient. 



Even at time, t= «> , a friction force acts on the plate — a situ- 
ation contrary to the previous case and the Rayleigh problem. 
Values computed by equation (22) are shown in figure 3. 
The force on the magnetic field is given by 


F{ unit area= 



(23) 


where u is given by equation (21). The force is finite in this 
case. 

Another case which may be of interest arises when the 
magnetic field is fixed on the plate but a pressure gradient 
compensates for the action of the magnetic field at y= <*> . A 
flow of this character may be imagined to exist in the early 
stages of boundary-layer formation on the walls of a two- 
dimensional channel. The conducting fluid is assumed to be 
pumped through a channel containing a stationary transverse 
magnetic field. A pressure gradient arises automatically to 
compensate for the resistance of the magnetic field. The 
equations then describe the boundary layer as a function of 
time, after the pump has started. Eor this case equation (7) 
reduces to 


5T+m, («-«.)+! g 


bhi 


The pressure gradient ^ ^ is found by examining the bound- 
ary condition mentioned above at y— CO . At y= CD 


Therefore, at y= <» J 


du dhi 

it~W~ u=0 


, 1 dp 

■ miU ~ + p*r° 


and the above equation becomes equation (12). Since the 
boundary conditions are the same as when the magnetic field 
was fixed relative to the fluid, the results shown in figures 
2 (b) and 3 apply here also. 

BOUNDARY LAYER ON A SEMI-INFINITE FLAT PLATE 

In the analysis of the boundaiy layer on a flat plate moving 
impulsively, it was found that the shape of the boundary- 
layer velocity profile changes with time, necessitating a new 
calculation at each instant. Therefore, a simple transforma- 
tion of the coordinates would not yield a similar solution 
wherein all profiles could be drawn as one. It is then to be 
expected that in the flow of an incompressible fluid over a 
semi-infinite flat plate in the presence of a magnetic field, a 
similar solution will not exist. Since the boundary-layer 
solution of Blasius is not known in closed form and the 
velocity profiles at two different stations along the plate are 
not similar, a numerical solution to a series expansion is 
found by the method used in reference 33. The first or 
zero-order term is found to be the same as the Blasius solu- 
tion which is well known and is tabulated, for example, on 
page 103 of reference 28. The first nonzero term will be 
calculated for the cases studied here. 

The boundary-layer velocity and temperature profiles will 
be found for three cases. The conductivity of the fluid is 
assumed to be uniform throughout the flow field in the first 
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and second cases. A particular variation is assumed in the 
third case. The magnetic field is fixed relative to the plate 
in the first and third cases and relative to the fluid in the 
second case. A fourth case is shown to be equivalent to the 
second. In all cases the fluid will be assumed to be incom- 
pressible and to have a Prandtl number of 1 .® 

MAGNETIC FIELD FIXED RELATIVE TO THE PLATE 

Velocity profile. — The magnetic field lines of force are 
assumed to be perpendicular to the free-stream direction and 
to begin at the leading edge of the plate as indicated in 
figure 4. An experimental setup may be imagined to consist 
of a wing or flat plate in a stream of conducting fluid such as 
mercury, sodium, or salt water. A magnetic field is impressed 
across the flow field using either a permanent or electro- 
magnet. 



Figure 4. — Fluid flowing past a semi-infinite flat plate in the presence 
of a transverse magnetic field. 


On the basis of the usual boundary-layer assumptions, 
equation (7) reduces to 

bu . bit . bhe ., 

U Tx +V b^ +miU=v b^ (24) 

The boundary conditions are that u=0 at y= 0 and u=u w 
upstream of the plate. At y= a >, bu/by=0 and bu/bx= — m t . 
Outside of the boundary layer the continuity equation 
requires that c )vlby=m l ; therefore, v—** as y— >». If the 
flow field is restricted to the vicinity of the surface of the 
flat plate, no real difficulty arises. In an actual problem 
the magnetic field will not extend uniformly to large distances 
from the plate. At the location where the magnetic field 
strength dies out another viscous layer will develop. This 
slip layer will be ignored. Equation (24) is the same as the 
incompressible boundary-layer equation except for the 
additional term mjTi. As was done in the Blasius boundary- 
layer solution, the transformation, 

(25) 

is introduced. The stream function, i p, is defined as 

f=-Ju~ni Uo^nfmxJi+mx^ 3 + (mx) 3f2 f t + (mx) 2 ^ . . .] (26) 

where the functions f 0 , f 3 , . . . are functions of 17 only. 

The series expansion in -\frixx is suggested by equation ( 21 ), 

* Recent theoretical computations by Hansen (reL 31) indicate that when the air is lonfaed 
up to 30 or 25 percent, a more realistic value for the Prandtl number would be about 0.3. 
When the gas becomes fully Ionised, the Prandtl number drops off to the order of 0.0L 



where the exact form of m is yet to be determined. If the 
expression ( 21 ) is expanded in a power series in mit, the odd 
power terms contributed by the first and second parts cancel 
each other. Also, if these terms are carried along in the 
following analysis, the corresponding functions are found to 
be zero. Therefore, only whole powers of mx will be 
considered in what follows. 

The velocity components and their derivatives are found 
from equations (25) and (26) as 

^b^b^ bif/by 
by £>£ by bn by 

u—u a [f 0 '+mxf 2 ' + {mxYfi' . . .] (27) 

b\f/_ bpbl; bpbrj 
v ~ bx d£ bx dij bx 

v=l^^fu;+mzf 1 '+(mx)\V+ . . .]- 

l^^fVo+Zmxj i +5(mx)'f<+ . . .] (28) 

Similarly, 

-l^[fo"+mxf 3 "+(mxyf i "+ . . .] + 

[2mr//+4(mr)V/+ . . .] (29a) 


lr tt °°VS V'" +mx f** + ( ma W+ • • •] ( 29b ) 

(mx)*fr+ • . •] (29c) 

When equations (27) through (29) are substituted into 
equation (24) and the like powers of 

<tB* M 2 

mx= — -x==- (30) 

are equated, the following differential equations result: 

(31) 

K'=lo , j> , -hoV-hdo n +1°' (32a) 


//'= 2/.'// (32b) 


The boundary conditions are 

jo=ii=j<= ... =0 at 7j=0 
//=//=//= ... =0 at y=0 
jo= 1 and//= — 1 at 7 /=“ 
//=/«'= . . . =0 at 71= m 
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The differential equations for / 2 and / 4 are linear ordinary 
differential equations with variable coefficients. Since the 
coefficients (Blasius solution) are not known in closed form, 
the equations must be integrated numerically. Only the 
functions/ 3 and/ 4 of the series will be found because the work 
involved is sizable and it is the intent of this study to find 
the over-all trends and gross effects rather than precise 
results. 

The Runge-Kutta method was used to integrate both the 
homogeneous and nonhomogeneous forms of equations (32a) 
and (32b). The integration is started at the surface of the 
plate using two boundary conditions and then assuming a 
value for /a" and//'. The correct solution is found by a 
proper combination of the numerical solutions to the homo- 
geneous and nonhomogeneous equations. The factor by 
which the homogeneous equation is to be multiplied is found 
from the boundary condition at t] =co ■ The numerical 
results for /a and/ 4 are tabulated in table I (a). The velocity 
profiles shown in figure 5 (a) indicate that the flow separation 
at mx=0.5 will probably not be predicted if a sufficient 
number of terms are taken in the series (27). The first 
three terms (//, //, and / 4 ') describe the velocity in the 
boundary layer quite well up to at least nuc=0.2. 

The skin-friction coefficient and the displacement 5* of an 
incoming streamline are 

_ 0.664-1.789 mz+0.706 m *x*- ... + ... , 

Cf 

8*=mxy+ (1.73+0.54 mr+1-34: . . . ) 

, (33b) 

The unit generating the magnetic field is assumed to be 
carried by the aircraft. The drag on this unit caused by 
the magnetic-field-fluid interaction is 

F/ unit area= J aBpu dy 

Substituting equation (27), 

, 

■F/unit area=crR 0 2 yu„vx {fo'+mxfi'-V • ■ ■ )dy 

Jo 







=aB 0 i ^u a vx (f a +mxf t + . . . ).”->” (33c) 

The force is finite only for a conducting layer of finite 
thickness. 

Temperature profile. — When the pressure gradients are 
reasoned to be everywhere zero and the usual two- 
dimensional boundary-layer assumptions are made, equa- 
tion (11) reduces to 


n b T, n bT t> a b*T. /buV , 0<n 

puC >bz +pv0p by <rBo ^'~Pr by i+ll \by) (34) 

If total energy is defined as 


E=C p T+^r 





(a) Velocity. 

(b) Temperature distribution for a cool wall, T„>= T a . 

(c) Temperature distribution for an insulated wall. 

Figure 5. — Boundary-layer profiles for fluid flowing over a semi- 
infinite flat plate with a transverse magnetic field fixed relative to 
the plate. 

equations (24) and (34) can be combined to yield a differential 
equation for the transport of total energy E; that is, if 
equation (24) is multiplied by u, added to equation (34), and 
the Prandtl number, Pr, is assumed to be 1, 

bE . bE b*E 

U b^ +V b^ =v W 


(35) 


(36) 
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TABLE I.— MAGNETIC FIELD FIXED RELATIVE TO PLATE 




(a) Stream function 



(b) Temperature func- 
tion for cold plate 

( 0 ) Temperature func- 
tion for Insulated plate 

V 

A 

A' 

A" 

A 

u 

A" 

Tt 

7Y 

u»V2(7, 

T 1 

Ua>*f2C r 

vy 

U<oV2 C, 

0 

0 

0 

-a 894 

0 

0 

0.353 

0 

-a 352 

0 

0 

.2 

-.018 

-.178 

—.888 

.007 

.069 

.330 

-.047 

— . 116 

.024 

.236 

.4 

-.071 

-.354 

-.888 

.027 

.132 

.287 

046 

.110 

.094 

.465 

.0 

-.159 

-.525 

, JCW 

.059 

.183 

.213 

-.001 

.330 

.209 

.679 

.8 

-.281 

-.687 

, 7K£) 

.099 

.216 

.117 

.087 

.538 

.304 

.867 

1.0 

-.433 

-.839 

-.730 

.144 

.228 

.068 

.213 

.715 

.554 

3.024 

L 2 

-.615 

-.978 

-.660 

.189 

.218 

-.111 

.370 

.859 

.771 

1.139 

1.4 

-.824 

-1. 102 

— . 579 

.230 

.184 

-.228 

.653 

.965 

1.000 

3.207 

1.6 

-1.055 

-1.209 

-.489 

.201 

.127 

-.338 

.754 

1.028 

1.250 

1.224 

1.8 

-1.300 

-1.298 

-.394 

.279 

.050 

-.431 

.962 

L047 

1.492 

1.187 

2.0 

— 1. 673 

-1.367 

-.296 

.280 

— . 044 

-.500 

1. 169 

1.020 

L722 

2.302 

2.2 

-1.852 

-1.416 

-.197 

.201 

-.148 

-.539 

1.367 

.954 

1.930 

.073 

2.4 

-2.138 

-1.446 

-.103 

.221 

-.257 

-.544 

L549 

.854 

Z 108 

.811 

2.6 

-2.429 

-1.458 

-.017 

.169 

—.363 

-.514 

1.707 

.730 

Z 253 

.628 

2.8 

-2.720 

-1.454 

.059 

.070 

—.460 

-.450 

1.840 

.593 

Z 359 

.440 

3.0 

-3.010 

-1.435 

.121 

-.024 

-.641 

-.358 

1.944 

.463 

Z 429 

.259 

3.2 

-3.294 

— L 403 

.168 

-. 139 

-.602 

—.247 

2.021 

.321 

Z4G4 

.096 

3.4 

-3.572 

—1.369 

.201 

-.204 

-.040 

-.125 

2.074 

.205 

Z469 

-.039 

3.6 

-a 841 

-1.327 

.218 

-.393 

-.653 

-.005 

2.105 

.108 

Z451 

-.143 

3.8 

-4.102 

— L 283 

.222 

-.623 

-.642 

.100 

2.118 

.033 

Z414 

-.213 

4.0 

-4.354 

-1.239 

.216 

-.649 

-. 611 

.199 

2.110 

-.020 

Z307 

-.253 

4.2 

-4. 598 

— L 197 

.200 

-.707 

-.564 

.269 

2. 112 

-.054 

Z315 

-.207 

4.4 

-4.833 

-1. 159 

.180 

-.874 

-.605 

.-314 

2.099 

-.072 

Z 262 

-2C0 

4.6 

-5.062 

— L 125 

.156 

-.968 

-.440 

.335 

2.083 

-.078 

Z 212 

-.239 

4.8 

-5.284 

-1.097 

.131 

-1.049 

-.373 

.334 

2.068 

-.075 

Z 107 

-.210' 

5.0 

-5.501 

-1.073 

.107 

—1. 117 

-.307 

.316 - 

2.054 

-.068 

2.128 

177 

5.2 

-5.713 

-1.054 

.085 

-1. 173 

—.247 

.285 

2.041 

-.058 

ZO90 

144 

5.4 

-5.923 

— L 039 

.065 

— L 217 

-.194 

.248 

2.030 

-.047 

Z070 

-.113 

5.0 

-0.129 

-1.028 

.049 

—1.251 

-. 148 

.208 

2.022 

-.037 

Z 050 

-.087 

5.8 

-6-334 

-1.019 

.036 

-1.277 

-.111 

.168 

2.010 

-.028 

Z035 

-.004 

6.0 

-0.537 

-1.013 

.026 

-1.296 

-.081 

.132 

2.011 

-.021 

Z024 

-.047 

6.2 

-6-739 

-1.009 

.018 

-1.309 

-.057 

.101 

2.007 

-.015 

Z016 

-.033 

6.4 

—0.940 

— L 006 

.012 

-1.319 

-.040 

.075 

2.005 

-.010 

Z011 

-.022 

6.0 

-7. 141 

-1.0Q4 

.008 

-1.320 

-.027 

.054 

.2.003 

-.007 

Z 007 

-.016 

6.8 

—7.342 

-1.002 

.005 

-1.330 

-.018 

.038 

2.002 

-.004 

ZOOS 

-.010 

7.0 

-7.542 

-L001 

.003 

-1.333 

—.012 

.020 

2.001 

.003 

Z 003 

-.006 

7.2 

-7.743 

-1.001 

.002 

-1.335 

-.008 

.017 

2.001 

.002 

Z 002 

-.004 

7.4 

-7.943 

-1.001 

.001 

— 1.33G 

-.005 

.011 

2.000 

.001 

Z 002 

-.002 

7.6 

-8.143 

-L0 

.001 

— 1.337 

-.003 

.007 

2.000 

.001 

Z001 

-.001 

7.8 

-a 343 

-1.0 

0 

-1.337 

-.002 

.004 

2.000 

0 

ZOOO 

-.001 

8.0 

-a 643 

-1.0 

0 

-1.338 

-.001 

.002 

zooo 

0 

ZOOO 

0 

8.2 

-a 743 

-L0 

0 

-1.338 

0 

.001 

2.000 

0 

ZOOO 

0 

8.4 

-8.943 

-1.0 

0 

-1.338 

0 

0 

ZOOO 

0 

2000 

0 

8.0 

-9.143 

-1.0 

0 

-1.338 

0 

0 

ZOOO 

0 

zooo 

0 


Although the temperature distribution is altered by the 
presence of the magnetic field because u and v are affected, 
the total energy of the conducting fluid is not. The kinetic 
energy removed by the force of the magnetic field is exactly 
equal to the heat generated by the electric current, inde- 
pendent of the field strength B 0 and the conductivity <r. An 
analogy may be drawn by considering a similar situation for 
an electric motor. Assume the motor to be rotating at a 
given speed. The power is then turned off and the arma- 
ture is short-circuited. The motor will decelerate at a given 
rate depending on the field strength, resistance of the circuit, 
and rotational speed, just as in the fluid flow problem. The 
temperature rise of the wires, insulation, etc., however, de- 
pends only on the initial kinetic energy of the rotor. If the 
conductivity of the wires (or fluid) is reduced, the time re- 
quired to stop the motion will be increased but the total 
f/cr or joule heat is not changed. 

A first-order estimate of the influence of the magnetic 
field on the temperature profile can be found from equation 
(34) together with equations (25) through (29) and table 
I (a). Assume that 

T=T 0 (r 1 ) + T,( v )mi + T i (v)nvh?+ . . . (37) 

where, T 2 , T t , . . . are functions only of y=y and 

Ti=T 3 =T 5 . . . =0 as in the part of the problem dealing 
with the velocity profile. The derivatives of the tempera- 


ture with respect to x and y are found as 


bT_bTblj . bT by 
bx btj bx~^by bx 

^=ir mxT i +2m*x?T i + . . . - 
b x x[_ 

l(T 0 '+mxT i '+TrMT i '+ . . . )] 

sf=51 §-■ < r "'+ ’’“’JV-WIY+ 


y ( 38 ) 


.) 


When the various expressions are inserted into equation (34) 
and terms containing like powers of mx are equated, the 
differential equations for T 0 and T 2 are 


To" , T e 'f 0 , uJUJ'Y 
Pr ‘ 2 ‘ C v 


(39) 


T" , Ti'fo 
Pr‘ 2 


fo'T i +±f i T 0 '+'^- (f,'*+2f o "f 2 ")=0 , (40) 


The boundary conditions are 


T 0 =T W at ij=0 

T 0 =T a at 77 = co 

T 2 "=r 2 '=0 at 7 ;= co 
T t = 0 at 17=0 
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The solution to equation (39) was first found by Pohlhausen 
and is discussed on page 246 ff. of reference 28. The func- 
tion T 0 (ij) is the temperature profile for a boundary layer 
on a flat plate and is given by (assuming Pr= 1)' 

(l-//)+f^ (l-// 1 ) (41) 

The derivative of T 0 to be inserted into equation (40) 
becomes 

To' = -So" (r w -T a So') (42) 


Equation (40) then becomes, for the insulated plate case 
(Pr=l), . * 



S.'Tr - ^ [| fJSfo"- (/ o , ) 1 ~2/ o ' / y/ / ]=0 

(46c) 


The temperature profiles (table I (c)) for several values of 
mx are shown in figure 5 (c). It is to be noted that the re- 
covery temperature is not altered by the magnetic field. 
Such a result is to be expected since the magnetic field and 
the fluid at the wall do not interact. 


For simplicity, the temperature of the plate will be assumed 
to be the same as that of the fluid far from the plate; that 
is, T w —T a . Equation (40) then becomes (Pr=l) 

T a "+§ TS-fST^ [| SJo" Q-/ o , )+(/o , ) 1 +2/ o ,/ / a ,, ]=0 

(43) 


Equation (43) was integrated numerically by the same 
method used for equation (32a) and the results are tabulated 
in table I (b). Several temperature profiles are shown in 
figure 6 (b) . The quantity of heat transferred to the plate , 
per unit time is 



(44a) 


The local convective heat-transfer rate is (Pr=l) 


h g __ p u - C p 
uJ/2C p %[W x 


(0.664-0.704 mx— ../+...) 

(44b) 


MAGNETIC FIELD FIXED RELATIVE TO THE FLUID 


Velocity profile. — In this case the magnetic lines of force 
move past the plate at the free-stream velocity. The fluid, 
which is decelerated by the viscous action of the plate, re- 
ceives a push from the magnetic field which counteracts the 
viscous effect. Such a situation exists when a plate, vane, 
or wing moves through a stationary electrically conducting 
fluid and magnetic field. The boundary layer on a wing of 
an aircraft flying in the ionosphere over the magnetic poles 
of the earth could be considered to correspond to this 
case. Under these circumstances, equation (7) reduces to 


bu . bu 


uBo 1 . . &u 

°-{u— 

p by* 


(47) 


If the transformation' of coordinates and the stream function 
described by equations (25) and (26) are introduced into 
equation (47), a set of ordinary linear differential equations 
with variable coefficients is found for f 2 , j), . . . 

§A/o"+//-l (48a) 


In the limit as the conductivity or field strength vanishes 
(m-> 0), equation (44b) reduces to 


fi)m=. 0 — 


fo"u„pC P 


(45a) 


or by virtue of equation (33a) 


m-Q — 2 


(45b) 


SS"=2j 0 'j^jSjS-\s o jS'-\w , -\j i So ,, +jS (48b) 

The boundary conditions are 

f 0 =fi t=ft= ... =0 at 77=0 

jo'- =//=//= • • • =0 at 77=0 

j 0 '=l and/i'=//= . . . =0 at 77=00 


Equation (45b) is often identified as Reynolds analogy. 
A similar relation of equal simplicity was not found for the 
problems treated in this paper. 

The differential equation for flow over an insulated plate 
is obtained from equation (40) by introducing the proper 
expression for Tj. The solution of equation (39) for the 
insulated plate case is (see, e. g., ref. 28) 

To=T«+gr ( W/ s) (46a) 

and therefore, 

T 0 '=-^So'So" (46b) 


The integration for f a and fi was carried out by the same 
method described in the previous section and tho results are 
tabulated in table II (a). Several velocity profiles are shown 
in figure 6 (a). The difference between the first- and second- 
order curves is small at small values of mx. However, the 
separation of the 7nr=0.5 curves indicates that the series 
expression (27) does not converge rapidly. More than two 
terms in the series must then be used if the velocity profile 
is to be predicted accurately at values of mx higher than 
about 0.2. The local skin-friction coefficient is 

0.664+2.293 mx— 2.768 m¥+ , 

v m (49a) 
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(a) Velocity. 

(b) Temperature distribution for a cool wall, T m . 

(c) Temperature distribution for an insulated wall. 

Figure 6. — Boundary-layer profiles for fluid flowing over a semi- 
infinite flat plate with a transverse magnetic field fixed relative to 
the fluid. 


and the displacement thickness is 
5*=(1.73 — 2.21 mx+5.22 

i 

(49b) 

The disturbance caused by viscous action on the plate 
results in a force on the magnetic field, which is given by 



F/unit area= 



o-J5 0 *(i4„— ^ u)dy 


From equations (27), (66), and table II (a), 

F/unit ar ea= aB 0 ^u m vxj o (1 —f,' — mxft — . ■ . )(h 

= (1 .73-2.21 mx+5.22m*x 4 — 

... + .. .)cB 0 i -yfu^vx 
r=cu a> 8*B 0 i 

The force would appear as a reaction on the unit generating 
the magnetic field and not on the model which disturbs the 
fluid. 

Equation (47) is identical with the differential equation 
for tie case when the magnetic field is fixed with respect to 
the plate and a pressure gradient exists just strong enough 
to compensate for the electromagnetic effects at tj= °o . 
The boundary conditions and, hence, also the numerical 
results of table II (a) and figure 6 (a) can be applied to this 
problem. 

Temperature profile. — Since the magnetic field lines are 
now moving relative to the plate, the expression for the 
electric current in the fluid is 


j=<xB 0 (u~uJ 

The equation for the transport of thermal energy in the 
boundary layer then becomes 


„dT, n bT nC v b*T, 

puCr bx+ pv0p bi~ aBo Tr bf 


(bu V 

W 


(50) 


If equations (47) and (50) are combined to describe the 
transport of total energy in the boundary layer (Pr=l) 


bB , bE, , . b*E 

PU -v a (p—uJ<rBo*=p 


(61) 


where, as before, E— G„T+ — • It is seen that the presence of 

the magnetic field adds energy to the local element by doing 
work on the fluid retarded by viscosity. The magnitude of 
the work depends on the velocity, electrical conductivity, 
and field strength. 

When the transformation of coordinates, equation (25), 
and the series expansion, equation (37), are inserted into 
equation (50), a set of differential equations is obtained for 
T 0} T 2 , T 4 , . . . The relation for T„ is the same as equation 
(39) and that for T 2 is 

+%T/-fo'T 1 +^ [(l-//)*+2/ "/ 2 "]=0 (52a) 
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TABLE II.— MAGNETIC FIELD FIXED RELATIVE TO FLUID 




( a ) Stream function 



( b ) Temperature func - 
tion for cold plate 

( c ) Temperature func - 
tion for Insulated plat © 

V 

/» 

/>' 

h" 

u 

u 


Ti 

Ti' 

tt c oty2Cp 

Ti 

«« >V2C r 

Ti' 

Oa > V2C ; 

0 

0 

0 

L 147 

0 

0 

- 1.384 

0 

2.022 

3.228 

0 

.2 

.022 

.210 

.953 

-.028 

-.275 

- 1.362 

.337 

1.370 

3.163 

-.631 

.4 

.081 

.382 

.773 

-.110 

— . 542 

- 1.300 

.555 

.821 

2.985 

— L 121 

.6 

.172 

.520 

.606 

-.243 

-.793 

- 1.205 

.672 

.370 

2.723 

— L 477 

.8 

.287 

.626 

.463 

-.425 

— 1.022 

- 1.080 

.709 

.014 

2.403 

— 1.710 

1.0 

.420 

.702 

.312 

-.650 

- 1.224 

-.932 

.684 

-.253 

2.047 

- 1.830 

1.2 

.666 

.752 

.188 

-.913 

— 1.394 

-.765 

.613 

-.438 

1.677 

- 1.851 

1.4 

.710 

.777 

.072 

- 1.205 

- 1.529 

-.684 

.513 

— . 549 

1.312 

- 1.786 

1.6 

.876 

.781 

-.028 

- 1.622 

- 1.626 

-.394 

.398 

-.596 

.667 

- 1.652 

1.8 

L 031 

.767 

-.113 

- 1.854 

- 1.686 

-.202 

.279 

—.589 

.654 

— 1.460 

2.0 

1. 181 

.737 

-.184 

- 2.193 

- 1. 708 

-.014 

.165 

— . 641 

.383 

- L245 

2.2 

1.325 

.694 

-.241 

— 2. 634 

- 1.692 

.165 

.064 

-.464 

.158 

- 1.006 

2.4 

1.458 

. 642 

-.282 

— 2.868 

- 1.643 

.327 

-.019 

-.369 

-.019 

-.766 

2.6 

1.581 

.582 

-.310 

- 3. 189 

- 1.563 

.467 

-.083 

-.269 

-. 149 

-.639 

2.8 

1.691 

.519 

-.323 

- 3.492 

- L458 

.580 

-.127 

-.172 

-.237 

-.337 

3.0 

1.788 

.454 

-.324 

- 3.771 

- 1.333 

. 66 $ 

-.153 

-.087 

-.286 

-.166 

3.2 

L 873 

.390 

-.314 

- 4. 024 

- 1.195 

.713 

-.163 

-.016 

-.305 

-.030 

3.4 

1.944 

.329 

-.296 

- 4.248 

- 1.050 

.732 

-. 161 

.036 

-.301 

.069 

3.6 

2.004 

.272 

-.270 

- 4.444 

-.904 

.722 

149 

.072 

-.280 

.134 

3.8 

2.054 

.221 

-.240 

- 4.610 

-.763 

.687 

-. 133 

.092 

-.249 

.171 

4.0 

2.093 

.176 

-.208 

- 4.749 

-.630 

.632 

-. 114 

.099 

-.213 

.184 

4.2 

2.124 

.138 

-.176 

- 4.863 

-. 611 

.565 

—.094 

.097 

-177 

.ISO 

4.4 

2.149 

.105 

-. 145 

- 4.954 

-.405 

.490 

-.075 

. 0S9 

142 

.168 

4.0 

2.167 

.079 

-.117 

- 5.028 

-.315 

.413 

-.058 

.077 

-.111 

.144 

4.8 

2.181 

.058 

-.092 

- 6.081 

-.240 

.339 

—.044 

.064 

-.084 

.120 

6.0 

2.191 

.042 

-.071 

- 5. 123 

-.179 

.271 

-.033 

.051 

-.063 

.097 

6.2 

2.198 

.030 

-.053 

— 6. 154 

-.130 

.212 

-.024 

.040 

-.045 

.076 

6.4 

2.203 

.021 

-.039 

- 6. 176 

-.093 

.161 

-.017 

.030 

-.032 

.057 

5.6 

2.200 

.014 

—.028 

- 6.192 

-.066 

.120 

-.012 

.022 

-.022 

.042 

6.8 

2.209 

.009 

-.020 

- 6.203 

—.045 

.087 

-.008 

.016 

-.016 

.030 

6.0 

2.210 

.006 

-.013 

- 6.210 

-.030 

.061 

-.005 

.011 

-.010 

.021 

6.2 

2. 211 

. 004 

-.009 

- 6.215 

-.020 

.043 

-.003 

.007 

-.008 

.014 

6.4 

2.212 

.002 

-.006 

- 6.218 

-.013 

.029 

-.002 

.005 

-.004 

.010 

6.6 

2.212 

.001 

-.004 

- 5.220 

-.008 

.019 

-.002 

.003 

-.003 

.006 

6.8 

2.212 

.001 

-.002 

- 6.222 

-.006 

.012 

—.001 

.002 

-.002 

.004 

7.0 

2.213 

0 

-.002 

- 5.522 

-.003 

.008 

0 

.001 

-. 001 

.002 

7.2 

2.213 

0 

1:88} 

- 5.223 

-.002 

.005 

0 

.001 

-.001 

.002 

7.4 

2.213 

0 

- 5.223 

-.001 

.003 

0 

0 

0 

.001 

7.6 

2.213 

0 

0 

- 6.223 

-.001 

.002 

0 

0 

0 

.001 

7.8 

X213 

0 

0 

- 5.223 

0 

.001 

0 

0 

0 

0 

8.0 

2.213 

0 

0 

- 5.223 

0 

.001 

0 

0 

0 

0 

8.2 

2 . 213 

0 

0 

- 5.224 

0 

0 

0 

0 

0 

0 

8.4 

2.213 

0 

0 

- 6.224 

0 

0 

0 

0 

0 

0 

8.0 

2.213 

0 

0 

- 5.224 

0 

0 

0 

0 

0 

0 


When the expression (42) is inserted into equation (52a) and 
it is assumed that T w — and Pr= 1, the result is 

37'+§ zy-//r a + 7^f| fJo" (§-//)+ 

' O p L4 J ( 52b ) 

(W/)H-2//7,"J=0 

The boundary conditions are 

T 2 =0 at 77=0 
Ti — 0 at 00 

The numerical integration of equation (52b) was carried 
out by the method used earlier and the results are tabulated 
in table II (b). The temperature profiles shown in figure 6 
(b) correspond to the same values of rax shown for the veloc- 
ity profiles in figure 6 (a). The temperature gradient and 
the heat transfer to the wall are, respectively, 

fL-^VS( o - 332 + 2 o22 ”“+-- •> 

< 0 - 664+4 - 0 « »*+••■> <«> 

The recovery temperature for an insulated plate is altered 
by the relative motion between the magnetic field and the 
fluid at the plate. The differential equation for the pertur- 


'bation profiles are obtained by inserting equation (46b) into 
equation (52a). 

•f Qni J - 1 / 2 

2V'+ f T 2 '-UT 2 -^~ / 2 ///o"+^r- [(l-/o , ) 2 +2//7/ , ]=0 

(54a) 

Several typical temperature profiles (table II (c)) are shown 
in figure 6 (c). The recovery temperature at the wall is 
changed due to the^/c or joule heating caused by the relative 
motion between the magnetic field and the fluid at the plate; 
that is, the recovery temperature at the wall is given by 

r r ),_o=T 0 +1.614^«wH- ....... 

Tr),-o=T„+^ (1+3.228 mx+ (54b) 

The heat transferred to the plate may also be written as 
2 =h(T r —T a ), where T r is the recovery temperature given 
by equation (54b) and T w is the wall temperature. Since 
the magnetic field increases T r as well as t, the heat-transfer 
rate is increased in two ways. 

MAGNETIC FIELD FIXED RELATIVE TO PLATE— FLUID OF VARIABLE 
CONDUCTIVITY 

Velocity profile. — Numerous problems in fluid dynamics 
exist wherein the electrical conductivity of the fluid or the 
strength of the magnetic field changes with distance from 
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the surface. The magnetic field strength can be varied 
somewhat by the design of the apparatus. It is possible to 
adjust not only the strength of the field but also the variation 
with distance from or along the surface within certain phys- 
ical limits. A first approximation to the shape and magni- 
tude of most magnetic fields would be a uniform magnetic 
field of constant magnitude throughout the flow field. Such 
a magnetic system will once again be chosen for this case. 

The electrical conductivity of the fluid would be expected 
to vary in a complicated manner in directions along and 
perpendicular to the surface and from one situation to the 
next- The constant value chosen for the conductivity in the 
two previous cases is not realistic. The conductivity will 
vary under actual conditions of flow in a boundary layer. 
This will surely change the skin friction and heat transfer 
from that found in the last two cases. A logical variation 
will be studied in this section to see how large an effect a 
variable conductivity will have on c f and h. 

It is to be expected that the conductivity will vanish at 
the outer edge of the boundary layer and reach a maximum 
somewhere near the surface. The equations which, can be 
used to determine the number of ions (and hence conductiv- 
ity) in a stream at high temperature or under circumstances 
wherein combustion is taking place, are cumbersome. 
Therefore, a simple, approximate, analytical relation for 
the electrical conductivity as a function of temperature 
between the plate and outer edge of the boundary layer 
was not found. The conductivity of air after the passage 
of a normal shock wave has been measured in a shock tube 
and reported on page 339 of reference 3. The conductivity 
is reproduced in figure 7. At Mach numbers above 20, the 
results are approximated by a straight line. 

In view of the lack of knowledge of the variation of <r in 
the boundary layer, a linear relation between the velocity 



Figure 7. — Electrical conductivity of air as a function of shock wave 
Mach number as found in reference 3. 


decrement and the conductivity will be chosen for this 
case. Therefore, it will be assumed that 



(55) 


where a 0 is the conductivity at the surface, where u=0. 
Whether equation (55) is realistic is uncertain. Inserting 
the expression (55) into equation (24) gives 


bu , bu . m,u . . bhi 

u^-+v—+^~ (u a -u)=vj- t 


bx by v K 




(56) 


As before, an ordinary differential equation for each of the 
functions f 2 , /<, . . . results when the variables in equations 
(25) and (26) are inserted into equation (66) and like powers 
of mx are equated. 

// ' ' =////-§ Ufo ' -\foh' ' +fo' (1 -//) (57a) 


// ' '—‘Zfo'fi ^ foft ' ' +/,' - 2 /,'/.' 

(57b) 

The boundary conditions are, 


jo=jt=j<= ... =0 when ij=0 
//=//=//= ... =0 when tj= 0 
//=l,/s'=//= ... =0 when rj=co 


The numerical integration of equation (57a) was carried 
out for/ 2 by the method used in the other cases and the 
results are tabulated in table III (a). Typical boundary- 
layer velocity profiles are shown in figure 8(a). The local 
skin-friction coefficient is 


0.664—0.627 mx±. ■ . . 
■yjEe x 


and the displacement thickness is 


S*= (1.73+0.97 mr+ 



(58a) 


(58b) 


The force on the unit generating the magnetic field is 
expressed by 

F=<rB 0 \lu a ,vxj^ (1 —fo'—mxfi'— . . .) (//+ 
'77ixf/+ . . . )d77/unit area 


The integral is evaluated through the use of equation (57a) 
and table m (a) as 

F=crB 0 ^u a vx (0.664+0.649 mx+ . . . )lb/unit area (69a) 
or, integrating from the leading edge to x, 

o 

-?■=- x<r o B 0 \fu^vx (0.664+0.390 mx-\- . . . )lb/unit width 

(59b) 
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(a) Volocity. 

(b) Temperature distribution for a cool wall, T w =*T m . 

(c) Temperature distribution for an insulated wall. 

Figure 8. — Boundary-layer profiles for fluid flowing over a semi- 
infinito flat plate with a transverse magnetic field fixed relative to 
the plate-fluid of variable conductivity. 


Equation (59a) or (59b) expresses a drag force on the plate 
which is to be added to the friction drag. In coefficient 
form, the local drag per unit area on the plate due to the 
magnetic field is 


(1.328+1.299 mx+ . . . ) (60a) 

Vfi’e* 

When equations (58a) and (60a) are added, the net local 
force coefficient to a first order in mx is 


c D =— L= (0.664+0.701 ffli+ . . . ) (60b) 

V#«r 


The interference of the magnetic field reduces the skin friction 
at the expense of a higher drag. 

Temperature profile. — The thermal-energy transport 
equation is reduced in the manner used in the previous 
cases. The electrical or joule heating is now described by 





Equation (11) is then 


pC,u^+pvC, by 


bT * 0 B* 




(u a —u)v}- 


= Pr 


b*T . 

by i+l1 



(61) 

\ 

If equation (56) is multiplied by u and added to equation 
(61), the resulting equation for the transport of total energy 
is identical with equation (36). Once again, the tempera- 
ture distribution is changed but the total energy of the 
fluid is not altered by the presence of a magnetic field. 

The change of variables if, 17 , the stream function, xfr, and 
temperature expansion given by equations (25), (26), and 
(37), respectively, are now introduced into equation (61). 
A set of linear ordinary differential equations for T 0 , T S) 
T«, . . . is obtained for this case when like powers of mx 
are equated. The expression for T 0 is once again the same 
as equation (39). The differential equation for T 2 is 


T a " , f 0 T/ 
Pr' 2 


-J/T a +p 3 T 0 '+^- [(W/)/.' 4 +2/"//1=0 

(62) 


If it is assumed that Pr= 1 and the expression for T 0 ' is 
inserted into equation (62), the result is 



■f/T»+ 


[§/,//' Q~//)+ (1 -//)//H-2//^"]-0 (63) 

Numerical integration of equation (63) was carried out 
by the method used for the previous cases. The results 
are presented in table III (b) and several temperature 
profiles are shown in figure 8 (b). The coefficient for the 
convective heat-transfer rate is given by 

(0.664—0.206 mx- ... + ...) (64) 


1 

2^IW X 


B20G97 — 00 S3 
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TABLE III— MAGNETIC FIELD FIXED RELATIVE TO PLATE— VARIABLE CONDUCTIVITY 


(a) First -order stream function 

(b) Temperature 
function for 
cold plate 

(o) Temperature 
function for 
insulated plate 

V 

/> 

a 

h" 

Tt 

iy 

Ti 

n 





u c „>(2C, 

tZco 5 /2CV 

Ma’PC, 

«oo s /2 C, 

0 

0 

0 

-0.313 

0 

-0.103 

0 

0 

.2 

-.006 

-.062 

-.307 

-.012 

-.021 

.008 

.082 

.4 

—.025 

-.122 

-.289 

-.000 

.056 

.032 

.158 

.6 

-.055 

-. 177 

-.262 

.009 

.122 

.071 

.221 

.8 

—.095 

-.226 

-.226 

.039 

.173 

.120 

.268 

LO 

. -. 146 

-.267 

-.185 

.077 

.207 

.176 

.295 

1.2 

-.202 

-.300 

-.139 

.121 

.223 

.236 

.300 

1.4 

—.264 

-.323 

-.092 

.165 

.220 

.295 

.283 

L6 

-.330 

-.337 

—.045 

.208 

.199 

.348 

.246 

1.8 

— .3$$ 

-.341 

0 

.244 

.103 

.392 

.193 

2.0 

-.466 

-.337 

.064 

.272 

.117 

.424 

.128 

2.2 

-.532 

-.325 

U77 

.290 

.064 

.443 

.066 

2.4 

-.596 

-.307 

.106 

.298 

.010 

.447 

-.016 

2.6 

-.656 

-.283 

.128 

.294 

-.041 

.437 

-.081 

2.8 

-.709 

-.256 

.143 

.282 

-.085 

.415 

-.138 

3.0 

-.757 

-.226 

.150 

.261 

-.118 

. 3S3 

-. 180 

3.2 

-.799 

-.196 

.150 

.235 

-. 140 

.344 

-.20$ 

3.4 

-.835 

167 

.144 

.206 

-.151 

.301 

-.221 

3.6 

-.866 

-. 139 

.134 

.175 

. -.162 

.256 

-.221 

3.8 

-.891 

-. 113 

.121 

.146 

— . 145 

.213 

-.210 

4.0 

— . 911 

-.090 

.106 

.118 

132 

.173 

-.191 

4.2 

-.927 

-.071 

.090 

.093 

— . 115 

.137 

167 

4.4 

-.940 

—.054 

.075 

.072 

-.097 

.106 

-.141 

4.6 

-. 949 

-.041 

.060 

.054 

-.079 

.080 

-.116 

4.8 

-.956 

-.030 

.047 

.040 

-.063 

.060 

-.092 

5.0 

-.961 

-.022 

.036 

.029 

-.048 

.043 

-.072 

5.2 

-.965 

-.015 

.027 

.020 

-.036 

.031 

-.064 

5.4 

—.968 

-.011 

.020 

.014 

-.027 

.021 

-.040 

5.6 

-.970 

-.007 

.014 

.009 

-.019 

.015 

-.029 

6.8 

-.971 

-.005 

.010 

.006 

-.013 

.010 

-.020 

6.0 

-.972 

—.003 

.007 

.004 

-.009 

.006 

-.014 

6.2 

-.972 

—.002 

.005 

.003 

-.006 

.004 

-.009 

0.4 

-.972 

-.001 

.003 

.002 

-.004 

.003 

-.006 

6.6 

973 

-.001 

.002 

.001 

-.003 

.002 

—.004 

6.8 

-.973 

0 

.001 

.001 

-.002 

.001 

-.002 

7.0 

-.973 

0 

.001 

0 

-.001 

.001 

-.001 

7.2 

973 

0 

0 

0 

-.001 

0 

-.001 

7.4 

-.973 

0 

0 

0 

0 

0 

0 

7.6 

-.973 

0 

0 

0 

0 

0 

0 

7.8 

-.973 

0 

0 

0 

0 

0 

0 

8.0 

-.973 

0 

0 

0 

0 

0 

0 

8.2 

-.973 

0 

0 

0 

0 

0 

0 


-.973 

0 

0 

0 

0 

0 

0 

8.6 

-.973 

0 

0 

0 

0 

0 

0 


The influence of the magnetic field on the heat-transfer 
rate and skin friction is not so large in this case as in either 
of the two previous cases (eqs. (33a), (44b), (49a), and 
(53)). From these results it is apparent that the change in 
the heat-transfer rate and the skin friction is strongly 
dependent on the manner in which the conductivity (or 
magnetic field strength) vary through the boundary layer. 
A magnetic field which changes in strength with distance 
from the surface and a variable conductivity are , to be 
expected in the real physical situation. 

The differential equation for flow over an insulated plate 
is obtained by substitution of equation (46b) into equation 
(62). 


mx will be dimensionless only if multiplied by a factor which 
relates the various dimensional quantities. When the cor- 
rect conversion factors are inserted, reference 27, 


m=1.529X10 


or 


m=6.432X10- 10 


...KIONM 
(£H(S) 


per inch (05a) 


per inch (65b) 


T 3 "+ f ~ -f.'Tr-^h Mo tf/'+TT Ki-/oO// 2 +2//7"]=o 

The temperature profiles for several values of mx are shown 
in figure 8 (c) and tabulated in table HI (c). As would be 
expected, the recovery temperature at the wall is not altered 
by the magnetic field because there is no relative motion 
between the fluid at the wall and the magnetic field. 

MAGNETIC PARAMETER 

UNITS 

The importance of the magneto-hydrodynamic effect 
depends on the size of the magnetic parameter m=<rB 0 i lpu at . 
The usual engineering units are: cr in mhos/inch; B 0 in lines/ 
square inch; p in slugs/cubic foot; and u m in feet/second. If 
each of these quantities is expressed as such, the parameter 


where p 0 = 0.002378 slugs/ft 3 is the air density at sea level. 

Another method of finding the conversion factor for m is 
to combine the electric motor equations (ref. 35) 


.F=8.8507X10~ 8 B c li lb 

(66a) 

E=B„luX 10 -8 volts=^ 

(66b) 

5^=8.8607X10- 1!I crJB 0 a 5-52? 
ul z in. 4 


force/7 3 (8.8507) (172S)10~ 10 oB? 

per inch 

Po(p/po )« 2 0.002378 (p/po)u 

aB * 

:6.432X10- 10 r f - per inch 

(p/Po)u m * 
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In words, the parameter m expresses the ratio of electro- 
magnetic force density to dynamic pressure. 

The value of m can easily vary over many orders of mag- 
nitude. However, the density p and velocity u associated 
with a particular vehicle are generally fixed within relatively 
narrow limits. It remains, then, to adjust the conductivity 
and magnetic field strength so that the desired characteristics 
are achieved. The limits for the magnetic field strength he 
between zero and a maximum (discounting astrophysical 
possibilities) of about lOfiines/sq in. The electrical con- 
ductivity depends on the fluid considered and can vary over 
wide limits from very small values for insulators to approxi- 
mately 10 fj mhos/inch for metals. The results of several 
investigations of the conductivity of gases will be discussed 
in the next section. 

ELECTRICAL CONDUCTIVITY 

The information on the electrical conductivity 6 of gases 
under the special conditions of this report is not large. 
However, it is known that the air around a missile will be 
a conductor if: (1) the energy imparted to the air is great 
enough to ionize it; (2) metal is ablated from the missile; 
(3) the missile is in the ionosphere where the air is ionized 
by radiation from the sun; or (4) other special devices are 
used to ionize the air, such as irradiation by liigh energy 
particles. 

The electrical conductivity of air behind a normal shock 
wave has been measured in a shock tube and is reported on 
page 339 of reference 3. The results, replotted in figure 
7, indicate conductivities of the order of 1 to 5 mhos/inch 
for air behind strong shock waves. 

Firings of ultraspeed pellets, references 36 through 39, in- 
dicate that another phenomenon will enter the problem of 
conductivity. Photographs and spectrographs taken with 
special cameras show that metal is eroded or ablated from 
the pellet and swept off the model into the wake where it 
burns. Since the metal is in a liquid spray or vapor form 
while it is in the vicinity of the pellet, the conductivity there 
would bo increased by a sizable amount. It would be diffi- 
cult to estimate the orders of magnitude to be expected. 
The tests show that the air can be made a conductor by 
ionization of its molecules and atoms and/or by metal vapor 
or spray eroded from the surface of the vehicle. It is pointed 
out in reference 40 that the latter effect predominates when 
meteors enter the earth’s atmosphere. 

The ah 1 in the ionosphere is partially ionized by radiation 
from the sun so that it is a conductor without addition of 
energy from the aircraft. The amount of ionization at 
various altitudes and the total number of particles are listed 
in tables 2 and 9 of reference 41. A rough estimate of the 
conductivity is found from reference 7 as 

mhos/inch (67) 

c Good electrical conductors arc always good heat conductors. The ratio may he expressed 
approximately as (ref. 7), 


where ft Is the coefficient of thermal conductivity, * Is Boltzmann's constant, and e Is the 
charge on an electron 


where e, l«, m„ and U, are the charge, mean free path length, 
mass, and velocity, respectively, of the electron and N 0 and N 
are the particle density at sea level and the altitude in ques- 
tion. The factor n is the electron density. The values for 
<r range from 10 -4 mhos/in. at 62 miles altitude to about 32 
mhos/in. at 250 miles altitude. 

The large variety in fuels makes it difficult to bracket the 
magnitude of the conductivity of combustion products but 
information is available for several cases. It was reported 
in reference 5 that the conductive of the jet from a 50- 
pound thrust acid-aniline rocket motor was measured as 
being about 10 -e mhos/in. It was added, however, that the 
addition of a small amount of an impurity such as sodium 
would strikingly increase the electron density. 

A discussion is given in reference 6 of similar experiments 
and theory on the combustion of hjulrocarbon fuels. The 
experimental evidence required that the number of electrons 
in the combustion products (10 11 to 10 14 per cc) be several 
orders of magnitude greater than predicted theoretically 
(10 s to 10 9 per cc). The maximum conductivity can be 
estimated by equation (67) to be roughly 10 -2 mhos/inch. 

Although the conductivities discussed are low compared 
with those of metals, they appear to be large enough so that 
appreciable magneto-hydrodynamic forces can be exerted. 
If the conductivity is too low for a specific case, it may be 
possible that alkaline salts or some other additive can be 
used to increase the number of free electrons and ions. 

APPLICATIONS 

The possibility of using a magnetic field to prevent a hot 
fluid from coming in contact with a wall which it would dis- 
solve or destroy was discussed in the introduction. The 
extent to which this can be accomplished is measured directly 
by the parameter TO=(cr^ p 2 /pw aj ). The larger the parameter 
m is, the more effectively the “fluid suspension” process can 
be carried out. The influence of the magnetic field on the 
fluid is expressed in equations (33), (44), (49), (53), (58), . . . 
by the product mx. The distance x represents an “adjust- 
ment distance” required before the fluid responds a given 
amount or adjusts to the magnetic lines of force. If the 
magnetic parameter is very large, the adjustment distance 
will be very small and the fluid will follow the magnetic 
lines of force quite closely as indicated in figure 9. The 
magnetic force restrains any motion across the lines of force 
because the inertia forces are insignificant in comparison. In 
this limit the magnetic fines of force are the same as invisible 
solid boundaries. 


u<o Streamlines 



Wires'-'' --Magnetic lines 

Figure 9. — Sketch of fluid of high electrical conductivity flowing'over 

a surface. 
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If the parameter m is very small the adjustment distance 
x becomes very large and the influence of the magnetic field 
becomes negligible. 

The numbers quoted in the last section for the conductivity 
and field strength indicate that most of the applications will 
lie between the two extreme limits. 

FLOW OVER A HIGH-SPEED VEHICLE 

Magneto-aerodynamic effects may be used to alter the 
external flow over aircraft. Examples which illustrate the 
changes in the boundary layer will now be discussed. It will 
be pointed out that the skin friction on an aircraft might be 
reduced by judicious use of a magnetic field either to inhibit 
transition to turbulent flow or to alter the boundary-layer 
profile. 

It has been shown in reference 13 that the Reynolds num- 
ber of transition is raised when a coplanar magnetic field is 
applied. The results of reference 13 can be used as a rough 
estimate for flows of the boundary-layer type when the par- 
ticular velocity distribution may be approximated by a para- 
bolic profile. This is permissible since a coplanar magnetic 
field does not change the velocity profile, which is not so for 
the cases considered in this paper. 

The results of reference 14 for a transverse magnetic field 
apply to a two-dimensional channel and should not be used 
for flows of the boundary-layer type because of the change 
induced by the magnetic field on the velocity profile. For 
the flow in the two-dimensional channel the change in the 
velocity profile was largely responsible for the resistance to 
turbulence. The same will not. always be true for other types 
of flow fields. It will then in general be necessaiy to carry out 
a stability calculation in order to determine whether a trans- 
verse magnetic field tends to stabilize or destabilize the 
boundary layer. A situation wherein the magnetic field 
destabilizes the flow field is presented in reference 26, 
warning one that the presence of a magnetic field does not 
always decrease the likelihood of transition. 

The numerical results of table I and the equations (33) 
and (44) may he used to estimate the effects of a magnetic 
field attached to an aircraft on boundary layers of uniform 
electrical conductivity. The model analyzed theoretically 
will probably never be duplicated in a physical situation. 
However, circumstances may arise wherein the conductivity 
and magnetic field are such that they vary only slightly or 
are almost uniform in a layer of finite depth. This might 
occur, for example, (1) in the boundary layer of high-speed 
aircraft; (2) in the region between the surface of a blunt 
body and its shock wave; (3) in local regions near an aircraft 
when material is being ablated from the surface; or (4) in 
the ionosphere wherein air is ionized by radiation from the 
sun (ref. 41) and is then a uniform conductor, provided the 
influence of the aircraft is small. The fact that the layer is 
of finite thickness causes the equations to be in error at the 
interface. However, if the conducting layer is relatively 
thick the equations would serve as a good estimate in the 
vicinity of the surface. 

The second semi-infinite flat plate case which was analyzed, 
equations (49) and (53), could be thought to correspond to 
flight over the magnetic poles of the earth or through some 


other magnetic field. The electrically conducting fluid 
could come about by the same processes mentioned in the 
previous paragraph. As before, the magnetic field strength 
and electrical conductivity will probably not correspond 
exactly to the model analyzed. A close approximation is 
obtained when an aircraft flies through the ionosphore over 
one of the magnetic poles of the earth (provided the influence 
of the aircraft on the electrical conductivity of the air and 
on the magnetic field is small). Since the earth’s magnetic 
field is small (R 0 « 4 lines/sq in. from ref. 42), effects of this 
type are generally negligible in the lower reaches of the 
ionosphere. At the higher altitudes the magnetic parameter 
is large enough to cause a change in the heat-transfer rate 
but at those altitudes the mean free path is so largo that the 
continuum flow analysis of this report is no longer applicable. 

If the air in the boundary layer is ionized by frictional 
heating, the conductivity of the air will change with distance 
from the surface (fig. 10). The actual variation to bo 
expected would be difficult to predict theoretically. The 



Figure 10. — Airoraft moving through the atmosphere at high spood. 

solution for a boundary layer •with variabl6 electrical con- 
ductivity found in a previous section of this paper (eqs. 
(58), (59), (60), and (64)) and tabulated in table HI can bo 
used as an approximate answer. Equation (58) shows that 
the s kin friction con be reduced by a sizable amount at the 
expense of increased total drag (see eq. (60)). The advan- 
tages of using a magnetic field lie in the possibility of reducing 
the aerodynamic heating, as expressed by equation (64). 

The preceding discussion dealt ■with the advantages that 
may be achieved by reducing the skin friction, heat transfer, 
or turbulence. Favorable effects may also be obtainod by 
increasing the displacement thickness of the boundary layer. 
Such a device, figure 11, could be used on the lower surface 
of a wing or body to induce forces for lift, control purposes, 
or drag. The magnitude of these forces is easily changed 
by increasing or decreasing the field strength which in turn 
increases or decreases the displacement thickness. 

Several additional problems arise when magneto-hydro- 
dynamic effects are used. The system generating the 
magnetic field will need to be cooled since its parts will 
dissipate power. Also, the weight of the unit required to 
generate- a magnetic field of sufficient strength may be too 
large to consider, or of such a size that any advantages aro 
nullified. If the slow moving fluid is hot enough to radiate 
energy, a sizable quantity of heat may be transferred oven 
though the hot fluid is not in contact with the surface. An 
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Figure 11. — Schematic diagram indicating use of magnetic field to 

build up the boundary layer on a wing for lift or control purposes. 

evaluation of this effect would require a more complete 
analysis than the solutions developed in this paper. 

Whether the applications just discussed can actually be 
made depends on the magnitude of the quantities entering 
the parameter m. Consider, for example, an aircraft moving 
through the upper atmosphere at 25,000 feet per second. 
From figure 7, a conductivity of not over 5 mhos/inch is to be 
expected. If p/p 0 is taken as 10 — 4 , and the above values are 
inserted into equation (65b), 

tiji 

m=6.43X10- 10 ^-^lO- 0 ^ 2 

In figures 5, 6, and 8, it is seen that the product mx should be 
of the order of 0.1 for a notable magneto-aerodynamic 
effect to be brought about. If a: is taken as about 100 
inches, the magnetic field strength, B„, should be around 
1000 lines per square inch so that the product mx is equal 
to 0.1. A magnetic field of such a strength is not unreason- 
able. 

FLOW IN PROPULSION UNITS 

The discussion on the reduction of skin friction and heat- 
transfer rate to the walls of an aircraft by electromagnetic 
forces applies as well to the internal flow of reaction engines. 
It may be added that a uniform magnetic field across the 
jet stream corresponds to the channel flow problems studied 
in references 8 to 12 and 14. In general, imposing such a 
field increases the skin friction, equation (49a), and lowers 
the thrust because of the interference of the magnetic field. 
A compensating influence is realized through a reduction in 
turbulence in the stream (see, e. g., refs. 9 and 14) which 
lowers the skin friction, heat transfer, and probably the 
noise level of jet engines. Shallow magnetic fields appear to 
be the most promising for reducing the heat transfer if the 
flow is nearly laminar. 

The cases and illustrations considered so far have assumed 
that an electric potential E has not been applied across 
the fluid. To do so would add another force term (see eq. 
(2)) to equations (7) and (11). If a source of power is 
available, a combined electric and magnetic field can be 
used to maneuver the fluid. A possible application would 
bo to use the combined fields to induce an electric motor 
action to obtain very high exhaust velocities. 


CONCLUDING REMARKS 

The laminar boundary-layer solutions presented in this 
paper for the flow of an incompressible electrically conducting 
fluid over a fiat plate indicate the changes that will be 
brought about by a transverse magnetic field. It has been 
found that the skin friction and heat-transfer rate are re- 
duced if there is no relative motion between the plate and the 
magnetic field. The skin friction and heat-transfer rate in- 
creased in the case where relative motion was assumed. In 
all the cases studied the total drag was increased. 

The product of the magnetic parameter m and an adjust- 
ment distance x is the quantity which determines the effec- 
tiveness of the magnetic field to change the heat-transfer 
rate. Other studies have shown that a magnetic field will in 
certain cases inhibit transition to turbulent flow and thereby 
retain a lower heat-transfer rate. 

All the cases studied, except one, assume that the con- 
ductivity and magnetic field did not change with distance 
from the surface. A more realistic case of a variable conduc- 
tivity was considered in one example where a special variation 
was chosen. The reduction in heat-transfer rate and skin 
friction for a given value of mx was not as large as when the 
conductivity was constant. This points out the importance 
of considering the way in which the conductivity and mag- 
netic field strength change with distance from the surface. 

Aaies Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Mar. IS, 1957 
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